function alpha=realHilbert(n)
H = hilb(n+1);
s = zeros(n+1, 1);
for j = 0 : n
    s(j+1) = (-1)^j;
end
t = zeros(n+1, 1);
for j = 0 : n
    for i = 1 : j
        t(j+1) = t(j+1) + (-1)^(i+j) / i;
    end
end
r = s * 0.69315 + t;
alpha = gausslineqn(H, r, n+1);
end